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The properties of future singularities are investigated in the universe dominated by dark energy 
including the phantom- type fiuid. We classify the finite-time singularities into four classes and 
explicitly present the models which give rise to these singularities by assuming the form of the 
equation of state of dark energy. We show the existence of a stable fixed point with an equation of 
state w < —1 and numerically confirm that this is actually a late-time attractor in the phantom- 
dominated universe. We also construct a phantom dark energy scenario coupled to dark matter 
that reproduces singular behaviors of the Big Rip type for the energy density and the curvature of 
the universe. The effect of quantum corrections coming from conformal anomaly can be important 
when the curvature grows large, which typically moderates the finite-time singularities. 

PACS numbers: 98.70.Vc 



I. INTRODUCTION 

The increasing evidence from the observational data 
indicates that (linear) equation of state (EOS) parame- 
ter w lies in a narrow strip around w — —1 quite likely 
being below of this value Q]. The region where the EOS 
parameter w is less than —1 is typically referred as a 
phantom dark energy universe. This is caused by the 
fact that when the phantom EOS is constructed in terms 
of a scalar field the corresponding kinetic term is cho- 
sen to have a wrong sign (negative kinetic energy). Of 
course, this is not the only possibility: the phantom-like 
value for w may appear from Brans-Dicke (BD) scalar- 
tensor gravity, from non-standard (negative) potentials, 
from the non-minimal coupling of scalar Lagrangian with 
gravity or even usual matter may appear in phantom- 
like form. Recent works |^ |^ have been devoted to 
the study of phantom cosmologies produced by different 
models. 

The existence of the region with w < —1 (if such a 
phase in the universe evolution indeed occurs) opens up 
a number of fundamental questions. For instance, the 
entropy of such universe is negative (or the characteristic 
temperatures should be negative). The dominant energy 
condition (DEC) for phantom matter is violated, as a 
rule. The phantom dominated universe ends up with a 
finite-time future singularity called Big Rip or Cosmic 
Doomsday (see Refs. Ia,ll|)- This last property attracted 
much attention and brought the number of speculations 
up to the explicit calculation of the rest life-time of our 
universe! In its turn, such a study motivated the math- 
ematical investigation of singularities when the DEC or 
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the strong energy condition (SEC) is violated (see also 
Ref. 0). However, the Big Rip singularity is character- 
ized by the growth of the energy and curvature invariants 
with a divergent scale factor at Big Rip time. The energy 
scale may grow up to Planck one, giving rise to the sec- 
ond Quantum Gravity era. Eventually, quantum effects 
become important near the singularity where they may 
moderate or even prevent the singularity 0, Q . 

In such circumstances, it is fundamentally important 
to understand the properties of singularities in phantom 
dark energy universe, to classify them and to search for 
realistic ways to avoid the singularities. The present ar- 
ticle is devoted to the investigation of the above circle of 
problems. We study phenomenological models in which 
the pressure density p of dark energy is given in terms of 
the function of the energy density p, i.e., p = —p — f{p). 
When the function f{p) is zero, the EOS parameter w 
is equal to —1 (cosmological constant EOS). We classify 
the types of future singularities and show explicit exam- 
ples of phantom dark energy which realize all mentioned 
types of singularities. The evolution of such dark energy 
universe at late times (near to singularity) is investigated 
both analytically and numerically. 

The paper is organized as follows. The general system 
of dark energy coupled with dark matter is considered 
in Sec. II. The corresponding background equations are 
written in autonomous form which permits to analyze 
the fixed points and attractor solutions. We also study 
the stability of fixed points against perturbations and 
show the existence of a stable critical point with w < —1. 
Sec. Ill is devoted to the classification of finite-time fu- 
ture singularities (four different types) and their explicit 
realization in terms of the universe with the equation of 
state: p = — p — /(p). We present the existence of singu- 
larities of much less explicit type than the standard Big 
Rip where the scale factor is finite at Rip time (includ- 
ing the type given in Ref. 01 )■ We construct examples of 
phantom dark energy universes in which all four types of 
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singularities appear. The violation of strong and domi- 
nant energy conditions is investigated as well. We also 
present one example of the model which admits the tran- 
sition from w<— ltow>— 1. In Sec. IV the relation 
between the EOS function f{p) and the appearance of 
singularities is studied. Sec. V is devoted to the numer- 
ical study of attractor solutions in dark energy models 
with non-relativistic dark matter. The phase plane anal- 
ysis confirms the structure of singularities and the sta- 
bility around the fixed points. 

In Sec. VI we present coupled phantom/fluid dark mat- 
ter models for several specific coupling functions. These 
models generalize some of the examples of the previous 
sections. We explicitly show multiple scalar field models 
that exhibit the similar singular structure as in the cou- 
pled phantom scenario. In Sec. VII the role of quantum 
effects is investigated for the types of singularities classi- 
fied in Sec. III. In accordance with previous attempts in 
this direction |^ IQJ it is demonstrated that quantum 
effects can moderate the finite-time singularities. Sum- 
mary and some outlook are given in the final section. 



II. AUTONOMOUS SYSTEM 

Let us consider a system with two fluids: (i) dark 
energy with an equation of state: p — p{p), and (ii) 
a barotropic perfect fluid with an equation of state: 
P = WmP- We wish to study the case of a nonrelativistic 
dark matter {wm = 0), but we keep our discussion in 
general for the moment. In a spatially flat Friedmann- 
Lemaitre- Robertson- Walker (FLRW) metric with a scale 
factor a, the background equations are given by 



p + iH{p + p) = -Q, 

Pm + iH{prn + Prn) = Q, 
^2 

H = ^{P+P + Prn + Pr,i) 



(1) 

(2) 
(3) 



where H = a/a\s the Hubble rate, and = SttG with 
Newton's gravitational constant G. Here a dot denotes 
the derivative with respect to cosmic time t. We ac- 
counted for the interaction term Q between dark energy 
and the barotropic fluid. We also have the constraint 
equation for the Hubble rate: 



y(p + Pm)- 



(4) 



In analogy with scalar field dark energy models, we in- 
troduce the following "kinematic" and "potential" terms: 



PK = {p+p)/2, pv = {p~p)/2. 



(5) 



For scalar fields one has p = ^e^'^ + ^^(0) and p — 
\f-<j)^ — V{(j)) (a normal field corresponds to e = 1 and 
a phantom to e = —1). We shall introduce the following 
dimensionless quantities: 



k'^Pk 



y ^ 



k'^Pv 
3iJ2 



(6) 



Then the above background equations can be written in 
an autonomous form as an extension of Ref. jTll | : 



dx 
ON 



-i^+p'ip)] 
+3x \2x+ (1 



Q 



2HpK_ 
Wjn){l-x-y)] , (7) 



dN 



3x 



-y 



2Hpv' 

+3y [2x + {1 + w„,)(l - X - y)] 



1 dH 
HdN 



3X+-{1+Wrn){l 



y) 



together with the constraint 



2 

Kj P'ln 

3iJ2 



= l-x 



y, 



(8) 
(9) 

(10) 



where N = In a and p'(p) is a function of p defined as 
p'{p) = dp/dp. Since K^pm/i'^H'^) > 0, the parameter 
range of x and y is restricted to be 

x + y<l. (11) 

When Wm — and Q = 0, these are simplified as 

ix[x-y-p'{p)] , (12) 

-^l-p'{p)]x + iy[l + x~y\, (13) 



dx 

diV ~ 
dy_ _ 
dN ~ 
1 dH 

HdN 



y)- 



(14) 



Let us consider a situation in which p'{p) asymptotically 
approaches a constant, i.e., p'{p) — > w. In this case the 
equation of state for dark energy is 



P 

w = — 
P 



y 



y 



(15) 



In the phase plane in terms of x and y, the trajectory 
corresponding to the constant is a straight line, i.e.. 



1 



y ■ 



w 



(16) 



Setting dx/dN = and dy/dN = in Eqs. ^ and 
((T^ . one obtains the following fixed points: (i) {x,y) = 
(0,0), (ii) {x,y) = [{l + w)/2,{l-w)/2] and (iii) {x,y) = 
(0, 1). The fixed point (i) corresponds to the barotropic 
fiuid dominant solution (flm ^1)- The point (ii) is the 
dark energy dominant solution, that is, JIde = x+y — > 1. 
In the case (ii) we have 



1 diJ _ 3 

HdN ^ "2'"^ 



■ w) 



(17) 



which means that the Hubble rate increases for u; < — 1. 
The class (iii) corresponds to the potential-dominant so- 
lution with r^DE = X + y ^ 1. 

One can investigate the stability of the system by con- 
sidering small perturbations Sx and Sy around the fixed 
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point {xo,yo), i.e., x = Xq + 5x and y = yo + 5y. By 
Eqs. (|12|l and H13|l . we obtain the linearized equations 



— = 3{2xo - yo - w)(5a; - SxQfJy , 
dA* 



(18) 



—Sy = 3(-l + yo + w)Sx + 3(1 + Xq - 2yo)6y , (19) 
dA* 



This can be written by using a matrix A4 : 



d f dx\ , , f Sx 
dNiSy =^[sy 



(20) 



One can study the stabihty of critical points against 
perturbations by evaluating the eigenvalues of the matrix 
A4. For the class (i) one has X — —3w and A = 3, 
which means that the solution is an unstable node for 
K7 < and a saddle point for w > 0. For the class (ii) 
we obtain A — 3w and A = 3(w + 1). Therefore the 
solution is a stable node for w < — 1, a saddle point for 
— 1 < w < 0, and an unstable node for w > 0. Then only 
the stable fixed point is the class (ii) with w < — 1. The 
class (iii) corresponds to the eigenvalues: A = —3(1 + w) 
and A = —3. Since = — 1 in this case by Eq. (|15|l . 
one has A = 0, —3. Therefore the fixed point (iii) is 
marginally stable. 

The above argument corresponds to the one in which 
p'{p) asymptotically approaches a constant w. We can 
expect that this argument may be applied to the case 
with a dynamically changing p'{p) by following "instan- 
taneous" fixed points as in Refs. (12i,il3(|. We will check 
this behavior numerically in Sec. V. 



III. MODELS OF FUTURE SINGULARITIES 

In the present section we will consider the dark energy 
universe models which contain finite-time, future singu- 
larities. It is clear that depending on the content of the 
model such singularities may behave in different ways. 
That is why it is useful to classify the future singularities 
in the following way: 

• Type I ("Big Rip") : For t ^ tg, a^oo, p ^ oo 
and IpI — ^ oo 

• Type II ("sudden") : For t ^ tg, a as, p ^ Ps 

and IpI — > oo 

• Type III : For t ^ t^, a ^ as, p ^ oo and \p\ oo 

• Type IV : For t ^ ts, a ^ ag, p ^ 0, |p| ^ and 
higher derivatives of H diverge. 

Here ts, as and ps are constants with Os ^ 0. The type I 
is so-called the Big Rip singularity which emerges for 
the phantom- like equation of state: w < —1. The type 
II corresponds to the sudden future singularity in Q at 
which a and p are finite but p diverges. The type III 
appears for the model with p = —p — Ap" [T^ . which is 
different from the sudden future singularity in the sense 



that p diverges. This type of singularity has been discov- 
ered in the model of Ref. (Q] where the corresponding La- 
grangian model of a scalar field with potential has been 
constructed. The type IV is a new type of singularity 
which appears in the model described below. 

In this section we consider the dark energy EOS char- 
acterized by 



P - 



(21) 



where f{p) can be an arbitrary function in general. Note 
that such EOS maybe equivalent to bulk viscosity, see 
Refs. 15s|. The function f{p) oc p" with a constant a 



was pro pos ed in Ref. ||^ and was investigated in detail 
in Ref. [l|. ^ 



When Q = in Eq. (Q, the scale factor is 



given by 



a = ao exp 



f{p)J ' 



(22) 



where ag is a constant. Here it is assumed f{p) does not 
vanish for all values of p. When f{p) = everywhere, the 
standard ACDM cosmology is recovered, so the above 
choice of EOS is a sound way to study the deviations 
from such cosmology. 

In this section we do not implement the contribution 
of a barotropic fluid, i.e., p,n ~ 0, Pm — and Q = in 
Eqs. C3)-(@l). Combining Eqs. (HI and igj) with (EB, one 
finds 



t 



dp 



which is used later. Eqs. © and Q give 



(23) 



^ = -^{p + 3p)^^[2p + 3f{p)] . (24) 

For example, for the specific choice of EOS: f{p) — 
—2/3/3 -|- po sin{pi / p) , with constants po and pi, the 
universe iterates the transition between the acceleration 
(a > 0) and deceleration {a < 0). 



A. The model of transition from u)> ltoit)<— 1 

In this subsection we study the model in which EOS 
changes from w > —1 to w < —1. 

When w is a constant, one has a ^ ^2/3(i+m)) g^j-^j 

H - 2/3(1 -I- w)t for w > -1, and a ^ {ts - 
or H ^ —2/3(1 -I- w){ts — t) for w < —1. It is interesting 
to understand what kind of f{p) admits the transition 
between the region with w > — 1 and that with w < — 1. 
Let us consider the following simple model 



a{t) = ao 



t 



ts~t 



(25) 



Here n is a positive constant and < t < ts- The 
scale factor diverges with a finite time {t — > ts) as in 
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the Big Rip singularity. Therefore tg corresponds to the 
hfe time of the universe. When t <^ tg, a{t) evolves 
as t", which means that the effective EOS is given by 
w = —1 + 2/(3n) > — 1. On the other hand, when t ^ ts, 
it appears w = — 1 — 2/{3n) < —1. 

From Eq. (|25(l the Hubble rate is given by 



H = n 



1 1 



t ts-t 

Then using Eq. |0J, we find 

3nVl 1 



(26) 



(27) 



Hence, both H and p have minima aX t = ta/'2. with the 
values 



^min — , ; p: 



_ 48n2 



(28) 



By deleting t in p, one obtains 



(29) 



Here the plus sign in Eq. H29(l corresponds to the region 
t > ts/1 and the minus one to the region t < ts/2. Com- 
bining Eq. with Eq. (QJ, it follows 



,(ri.±|£ i-l2(^ 



3n 



(30) 



Therefore the EOS needs to be double-valued in order 
for the transition to occur between the region w < —1 
and the region w > —1. Such a double- valued EOS often 
appears when there is a first-order phase transition. As 
we will see later in Sec. VI, if there are two kinds of 
matter or energy, a value of p can correspond to several 
values of p. We should also note /(pmin) = 0, that is, the 
minima of p and H correspond to the transition point: 
w = —1. 

The singularity a,t t — ts corresponds to the Big Rip 
type characterized by a ^ oo, p ^ oo and \p\ — *■ oo for 
t ^ tg. In this region /(p) behaves as /(p) ~ 2p/3n due 
to Eq. (|30|) , which means that the pressure p is linear in 
p, i.e., p = —p— 2p/3n. Therefore this gives the constant 
EOS: w = p/p — —1 — 2/3n. In another asymptotic limit 
t 0, one also obtains the constant value of w, i.e., 
w = -1 + 2/3n. 

For the general case where w crosses —1, since w = 
— 1 corresponds to /(p) = 0, in order that the integral 
J dp//(p) in Eq. (|22|) is finite, /(p) should behave as 



f{p)-fo{p-por 



< s < 1 



(31) 



where the condition /(po) = is assumed. Since < s < 
1, /(p) should be multi- valued near p = po in general. 



B. Specific model 

In what follows we shall investigate a model character- 
ized by 



fip) = 



ABp' 



a+0 



Ap" + BpP 



(32) 



where A, B, a and /3 are constants. As is shown below, 
this dark energy scenario contains a rich structure from 
the viewpoint of singularities. 
If a is larger than /3, we have 



f{p) 



Ap" when p — s- 
Bp^ when p — )■ oo 



When a,/3 7^ 1, Eq. ^ gives 



a = ao exp 



{a - 1)A {f3 - 1)B 



(33) 



(34) 



For 1 > a > if A, i? > {A, B < 0), a has a minimum 
(maximum) oq at p = and a goes to infinity (vanishes) 
when p — > oo. For a>l>/3, ifyl<0 and B > {A > 
and B < 0), a has a minimum (maximum) at the non- 
trivial (non-vanishing) value of p and a goes to infinity 
(zero) when p vanishes or goes to a positive infinity. For 
a > 1 > A,B > {A,B < 0), a goes to infinity when 
p oo (p ^ 0) and a vanishes when p ^ (p ^ oo). 
For a > f3 > 1 , a goes to ao when p — > oo . Furthermore if 
A > {A < 0), a ^ {a oo) when p ^ 0. UA,B>0 
{A,B < 0), a is a monotonically increasing (decreasing) 
function of p. If A > and B < {A < and B > 0), 
a has a nontrivial maximum (minimum) at a finite value 
of p. 



C. The model ((sH) with a = 2/3 - 1 

In what follows we shall concentrate on the case that 
a = 2/3-1. (35) 
Then Eq. H34|l may be solved with respect to p: 



fl±(^-6A{P~l)ln^ 



(36) 



which is vahd for ^9 7^ 1. One has /(p) = ABp/{A + B) 
for /3 = 1, which means that the EOS becomes a usual 
linear equation p — wp with constant w, that is, w = 
— 1—AB/{A+B). As this case has been well investigated, 
we will not consider it in this paper. 

When a = 2/3 — 1, the pressure p is given by 



ABp'P- 
ApP-^ + 



B 



(37) 
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Now Eq. (|34|l can be rewritten as 

n-2(/3-l) 



a = Uq exp 



^(/3-l) 



(38) 



Eq. H38|l tells that if /3 > 1, a ^ oq when p ^ oo, and 
a ^ (a ^ oo) when p ^ and A > (A < 0). On 
the other hand, if < 1, a — )■ ao when p 0, and a — > 
(a — > oo) when p ~^ oo and yl < (A > 0). Furthermore 
Eq. ^ tells that if /3 > 1, 



-^-p-BpP , 



when p 
when p 



, (39) 
oo . (40) 



Therefore w = p/p ^ —1 — (—1 + 0) when p ^ and 
A > {A < 0) and w — s- +cx) (— oo) when p — > oo and 
B < (i? > 0). Thus, except for the case where p 
and A < 0, the dominant energy condition (DEC), which 
requires. 



p>0 , p±p>0 , 



(41) 



is violated. On the other hand, the strong energy condi- 
tion (SEC), 



p + 3p>0, p + p>0, 



(42) 



is not violated when p ^ oo and B < 0. We also note 
that if B < 0, the weak energy condition (WEC), 

P>0, P + P>0, (43) 

and the null energy condition (NEC), 

p + p>0, (44) 

are not violated, either. 

Moreover, if /3 < 1, one gets 

p —p — Ap'^^~^ , when p — > oo , (45) 
p —p — Bpf^ , when p ^ . (46) 

Therefore w +oo (— oo) when p and B < {B > 
0). Then DEC is violated but if i? < 0, the violation of 
SEC does not occur. When p oo, one has w ^ —1 — 
(-1 + 0) for A > (A < 0). Then DEC is not violated if 
A < although SEC could be violated. The violation of 
WEC and NEC does not occur in this case. 
From Eq. (|37() it follows that p diverges when 



(47) 



which occurs for A/B < 0. Then we have p +oo (— oo) 
ioT B > {B < 0) as p ^ Ps + 0, and p +oo (— oo) 
for S < (S > 0) as p ^ Ps - 0. Then when p^ ps + 
with B>0oip^ps~0 with B < 0, only the violation 
of DEC occurs, while other energy conditions are not 
violated. 



By integrating Eq. (|23|) . we find 



2 2A 

p 2 



4/3-3'' (2/3-1)5' 
= -VSnAit - to) = T, 



(48) 



where to is an integration constant. This is valid for 
/3 ^ 1, /3 7^ 3/4, and (3^1/2. When /3 = 3/4, instead of 
Eq. (gHl, we obtain 



In 1 r = T , 



Po Bp 



1 — ' J 



(49) 



where po is a constant introduced for a dimensional rea- 
son. For P — 1/2, it follows 



B Po 



(50) 



We use these relations together with Eqs. H36|) and H37f) 
in order to study the properties of singularities. 



D. Sudden future singularity for the model (1321 
with A/B <0 

Let us study the property of the singularity at p = 
Ps = {—A/B)^^/^^~^\ This corresponds to a finite value 
of the scale factor by Eq. (|^ . When p ~ Ps, Eqs. 
(|^ and (|5n|) give finite values of t, i.e.. 





4fJ-3 
2(3-1) 


4(/3- 


1) 






(4/3 -3) (2/3-1) ' 


for /3 ^ 3/4,1/2, and 








T« =-4 




Bpo)^'] 





for /3 = 3/4, and 



2A 

"b 



1 - In - 



A 
B^ 



(51) 



(52) 



(53) 



for /3 = 1/2. 

Since p is finite, the Friedmann equation Q also shows 
that H is finite but from Eq. H24|l . we find that a diverges 
because of the divergence of p and the scalar curvature 
R also diverges since R — 2k^ [p — 3p). From the above 
arguments the singularity described by Eq. H47|l corre- 
sponds to the sudden future singularity (type II) in 0. 



E. Classification of singularities for the model (I32II 

In the previous subsection it has been shown that for 
all values of /3 the sudden future (type II) singularity 
appears for the model (|32|l if A/B < 0. In addition to 
this there exists a wide variety of singularities depending 
on the values of /3. 
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1. p>l 

When (3 > 1 Eq. implies that p oo as r ^ +0 
(-0) for A/B > {A/B < 0). From Eqs. §^ and ^ 
we find a — > ag and |p| — s- oo as p ^ oo. Therefore this 
corresponds to the type III singularity at which p and \p\ 
diverge with finite t and a. The equation of state behaves 
a,s w —> —1 — Bp^^^ for p ^ cx) by Eq. (|4()|l . which means 
that w +00 (-oo) for B < (B > 0). 

Note that r c» as p ^ by Eq. For A/B > 

the value of r is restricted to be < r < +oo. When 
A/B < 0, T is restricted to be in the range t > Ts with 
negative [see Eq. (|51f) ]. From Eq. (|48|l we find dr/dp — 
at p — ps- Then two branches appear for the region of 
p, i.e., < p < ps corresponding to < t < +oo and 
Ps < p < +00 corresponding to < t < 0. 

2. 3/4 < /? < 1 

Let us next consider the case with 3/4 < /? < 1. We 
find p ^ oo as T ^ +0 from Eq. Eq. §^ tells 

that a ^ oo for A > 0, and Eq. (|45ll gives \p\ oo for 
p — !■ oo. This means that the singularity at t = to is the 
type I (Big Rip). Since w — > — 1 — Ap^^^"^^ for p ^ oo 
by Eq. we have w -1 + (-1 - 0) for A < 

{A>0). 

In the limit p ^ 0, one has r — > +oo (— oo) for A/B > 
{A/B < 0). The allowed region of r is < t < +oo for 
A/B > 0, whereas r < with positive for A/B < 0. 
In the latter case there exist two branches, < p < 
Ps corresponding to — oo < r < r, and ps < p < oo 
corresponding to < r < Ts. 

5. 1/2 < /3 < 3/4 

In the case of 1/2 < /3 < 3/4, Eq. gHl) iniphes that 
p ^ as r ^ +00 (-oo) for A/B > {A/B < 0). 
Meanwhile we have r— >— ooasp— s-oo. Then the value 
of T is not limited (— oo < r < +oo) for A/B > 0. We 
find — oo < T < r, for A/B < 0, in which case there exist 
two branches: < p < ps and ps < p < oo. Except 
for the type II singularity at r = for A/_B < 0, finite- 
time singularities do not exist for 1/2 < /3 < 3/4. This 
property also holds for /3 = 3/4 and (3 = 1/2, as can be 
checked by Eqs. and 

4. P<l/2 

When (3 < 1/2, one has p ^ oo as t ^ — oo and 
p ^ as T ^ -0 (+0) for A/B > {A/B < 0). Then 
the range of r is — oo < r < for A/B > 0, and < 
T < Ts corresponding to < p < p^ and — oo < t < 
corresponding to ps < p < oo for A/B < 0. 

In the Hmit p ^ or t — > 0, Eq. gSJ tells that p — > 
for < /? < 1/2, p ^ -B for f3 = 0, and p +oo (-oo) 



for p < and B < {B > 0). Therefore there exists a 
type II singularity for /3 < (a — > ao, p ^ 0, |p| — > oo as 
t ^ to). In this case we have w — + +oo (— oo) for _B < 
{B > 0) by Eq. 

In the case of < /3 < 1/2 one has Ina/oo oc t^~^^ . 
Since 1 - 1/(2/3 - 1) > 2 for < /3 < 1/2, and H are 
finite. However d^H/dt" diverge for n > —1/(2/3— 1) as 
long as 1 — 1/(2/3— 1) is not an integer. This corresponds 
to the type IV singularity in which higher derivatives 
of H exhibit divergence even if a, p and p are finite as 
t to. In this case w +oo (— oo) for i? < (B > 0) 
by Eq. gHl- 

When /3 = 0, we find a ~ aoexp[BT'^ /{UA^)] for 
p ^ 0, showing the absence of the singularity for higher 
derivatives of H. Therefore there is no any future singu- 
larity. In this case w — > -l-oo (— oo) for B < {B > 0) as 
p^O. 

The obtained results are summarized as follows: 

• For A/B < there is always the type II singularity 
irrespective of the values of /3. 

• Irrelevant to the sign of A/B, the types of singu- 
larities are different depending on the values of /3. 

1. I3>1: 

There is a type III future singularity. DEC is 
broken, w +oo (-oo) for B < {B > 0). 

2. 3/4 < /3 < 1: 

There is a type I future singularity for A > 0. 
DEC is broken for A > 0. w -> -1+0 (-1-0) 
for A<0 {A> 0). 

3. 1/2 < ^ < 3/4: 

There is no a finite future singularity. 

4. < /3 < 1/2: 

There is a type IV future singularity, w 
+00 (-oo) for B <0 {B>0). 

5. /3 = 0: 

There is no finite future singularity, but when 
p^O, w ^ +00 (-oo) for B < (B > 0). 

6. /3 < 0: 

There is a type II future singularity. DEC 
is broken but SEC is not broken for B < 0. 
w +00 (-oo) for B < (B > 0). 



IV. GENERAL STRUCTURE OF THE 
SINGULARITIES 

In this section we discuss the general structure of singu- 
larities. Especially, the relation between the singularities 
and the behavior of /(p) is clarified. 
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A. Type I and III singularities 

Let us investigate a situation in which p goes to infin- 
ity. Then the Hubble rate diverges from the Friedmann 
equation 10}, which leads to the divergence of all curva- 
tures. As in Eq. H3^^(l . we assume the form 



for p 



(54) 



Here B and /3 are constants. When (3=1, the EOS 
is a usual linear equation p = wp with a constant w(= 
-l-B). 

The model with f{p) = BpP was proposed in Ref. 6] 
and further investigated in Ref. [I4j . Since we are now 
interested in the structure of the singularities, we assume 
the behavior (|54|l with /3 > only when p goes to infinity. 
Then by using Eq. (EH), it follows 



a ~ ao exp 



,1-/3 



3(1-/3)5 



(55) 



When /? > 1, the scale factor remains finite even if p goes 
to infinity, which corresponds to the type III singularity 
[T3 |. When /? < 1, we find a — oo (a ^ 0) as p ^ oo. 
for 5 > (S < 0). 

Since the pressure is now given by 



p^-p-BpP, 



(56) 



p always diverges when p becomes infinite. If /3 > 1, the 
EOS parameter w = p/ p also goes to infinity, that is, 
w +00 (-oo) for B < (S > 0). When /3 < 1, we 
have w ^ -1 (-1 - 0) for S < (S > 0) as p ^ oo. 
By using Eq. (|23|l for the function (|54|l . one finds 

2 P-/3+1/2 1 

^-^0+ „ , , for /3^-, (57) 



and 



VSkB 1-2/3 

In p ^ 1 

^^*°+vf;;^' ^=2 



(58) 



Therefore if (3 < 1/2, p diverges in an infinite future or 
past. On the other hand, if /3 > 1/2, the divergence of p 
corresponds to a finite future or past. 

From the above argument, one can classify the singu- 
larities as follows: 

1. /3 > 1: 

There exists a type HI singularity. DEC is broken. 

w-^ +00 (-oo) if B <0 (B > 0). 

2. 1/2 < /? < 1: 

There is a type I future singularity for B > 0. DEC 
is broken in this case. When B < 0, since a ^ 
as p — > oo, if the singularity exists in past (future), 
we may call it Big Bang (Big Crunch) singularity, 
w ^ -1 (-1 - 0) if B < (B > 0). 



3. < /3< 1/2: 

There is no finite future singularity. 

When /3 < 0, it was shown in Ref. [Ti | that the type II 
singularity appears as p ^ 0. In the next subsection we 
shall investigate a more general model. 

B. Type II singularity 

Lets us consider the type II singularity at which p is 
singular but p is nonsingular at a finite time t = ts as in 
Ref. 0. The starting function /(p) is given by 



/(p) ^ C{po - p)- 



(59) 



where 7 is a positive constant. We concentrate on the 
case where p is smaller than po, but the situation is ba- 
sically the same for p > po by considering the function 
/(p) = C(p — po)~^- In the limit p — po, the pressure p 
becomes infinite because of the divergence of /(p). The 
scalar curvature R diverges since R = 2k^ (p — 3p) . The 
EOS is 



-1 - 



C 



(60) 



00 for C < 



which gives w —> —00 for C > and i 
as p ^ Po. 

From Eq. (|22|l the scale factor is integrated to give 



a ~ oq exp 



(po - pr+' 



3C(7 + 1) 



(61) 



which means that a is finite for p = po. Since the Hubble 
rate H is nonsingular by Eq. Q), d remains to be finite. 
On the other hand Eq. (|24|l implies that a diverges for 
p ^ Po- By using Eq. (|23|l we find the following relation 
around p ~ po: 



t^to 



{Po - p) 



7+1 



(62) 



where to is an integration constant. Then t is finite (t = 
to) even for p = po. The above discussion shows that the 
function /(p) in Eq. H59|) gives rise to the sudden future 
singularity. 

This type of singularity always appears when the de- 
nominator of /(p) vanishes at a finite value of p. The 
model (|54|l with negative /3 is a special case of the model 
(P|l with Po = 0. 



V. ATTRACTOR SOLUTIONS 

In this section we account for the contribution of a 
barotropic perfect fluid and investigate the attractor 
properties for several models presented in Sec. HI. We 
study the case of a non-relativistic dark matter {wm = 0) 
without a coupling between dark matter and dark energy. 
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A. Constant p'{p) with p'{p) 7^ — 1 



Let us first investigate the case in which p' (p) asymp- 
totically approaches a constant, p'{p) w, with w 7^ —1. 
One example is presented in Eq. (|30|l . i.e., 



p = -p± 



2p 

in 



1/2 



1/2 



(63) 



the 



where pc = 48n /(k t~). In the limit p ^ 00 
minus sign in Eq. (|()!^ gives the constant EOS, w = 
— 1 — 2/(3n) < —1, whereas the plus sign gives w — 
-l + 2/(3n) > -1. 

When — 1 < w < we showed in Sec. II that the fixed 
point {x,y) — [(l-t-w)/2, (l — u')/2] is a saddle point. One 
can easily find that this fixed point is the intersection of 
the line (|16|l and the line x + y = 1. For the model H63() 
with a plus sign, one has y/x + 1 — in/ \/l — (pc/ pY^^ 
by Eq. ((T3|) . Since \/l — [pc/ pY^"^ < 1, the range of x 
and y in the phase plane is restricted to be ?/ > (3n— l)a;, 
y — X and x>Q. Note that the line (I16II corresponds 
to y = (3n — l)x by setting w = —1 + 2/(3n). Therefore 
the intersection of two lines y = (3n — l)x and y = \ — x 
is a saddle point given by [x.y) — [l/(3n),l — l/(3n)]. 
This is plotted as a point A in Fig. ^ 

We run our numerical code for initial conditions of x 
and y that belong to the allowed range explained above. 
We find that the solutions gradually approach x — 
which gives the EOS: w = — 1 by Eq. l(T3|l . The point 
A corresponding to t« = — 1 + 2 /{in) is not an attractor 
point. As the trajectories approach a; = 0, the energy 
density p decreases to a minimum value pc- Figure 
shows that there is some portion of the energy density of 
a barotropic fluid (fim = l—x—y > 0) around a; = if the 
initial condition of x is much smaller than y. When the 
trajectories get closer to the line x + y = 1 in the middle 
of the way, the final attractor point is {x,y) — (0,1), 
corresponding to a universe described by a cosmological 
constant with a vanishing energy density of a barotropic 
fluid. 

When w < —1 there exists one stable fixed point: 
{x,y) = [(1 + w)/2, (1 — w)/2]. In scalar-field dark en- 
ergy models, the kinematic term is negative {x < 0) at 
this critical point. For the model l|()3|l with a minus sign, 
one can easily show that the allowed range of x and y 
is described by y > — (3n -I- l)x, y < 1 — x and x < 0. 
The intersection of the two lines y — —{in -I- l)x and 
y = 1 — X, i.e., {x,y) = [— l/(3n),l -t- l/(3n)], represents 
the above stable fixed point. In Fig. [21 we plot the phase 
plane in terms of x and y for the model H63() with a minus 
sign. The solutions approach the above stable fixed point 
at which the EOS is given hy w = —1 — 2/(3n) < —1. 
We recall that this attractor point corresponds to the Big 
Rip singularity with a divergent energy density at finite 
time tg. If we connect two functions (|63|l at the mini- 
mum value of p, the solutions starting from the region 
a; > approach the line x = 0, enter the region a; < 0, 



>~i 0.60 




0.1 

X 



FIG. 1: The phase plane for the plus sign of the model H63|l 
with n = 2. The allowed range corresponds to y < 1 — x, 
y > 5x and a; > 0. The point A is a saddle point given by 
X = l/(3n) — 1/6 and y = 1 — l/(3n) — 5/6. The solutions 
approach the region a:: = 0, corresponding to the equation of 
state: w = —1. 




-0.1 
X 



FIG. 2: The phase plane for the minus sign of the model H63|l 
with n = 2. The allowed range corresponds to y < 1 — x, 
y > —7x and a; < 0. The point B is a stable node given 
by a; = -l/(3n) = -1/6 and y = l-f l/(3n) = 7/6, which 
corresponds to the late-time attractor. The equation of state 
at this point is a constant: w — —1 — 2/(3n) — —4/3. 
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and then approach the above stable fixed point. It is in- 
teresting to note that the final attractor point is the Big 
Rip type even if the system begins from a non-phantom 
stage {x > 0). 



B. Ip'(p)I^cx) 

Let us next consider the case in which |p'(/o)| diverges 
asymptotically. One example is given by the model H37() 
with /3 > 1. In this case there exists a type III singularity 
at which p and p exhibit divergence with finite time and 
scale factor. In the limit p — > oo the pressure p is given by 
Eq. which means that p'{p) ~ —1 — f3BfP~^ —oo 

and w ~ — 1 — Bp^^^ —oo for B > 0. 

We shall study a situation in which the system starts 
out from the region x < (i.e., negative w) and then 
evolves toward w —oo. By Eqs. H15() and (|37|l one ob- 
tains y/x = ~l- 2{Au + B)/{ABu'^), where u = p'^^'^ > 
0. When A and B are both positive, we find y > —x for 
X < 0. Therefore the allowed range is restricted to be 
y > —X, y < 1 — X and a: < in the phase plane. 

In Fig. O we plot the trajectories of the solutions for 
the model ^ with /3 = 1.1, A = 2 and S = 1. We 
find that all solutions approach a; ^ — oo and y — > +oo 
along with the line x + y — 1. The points — 
[{l+p'{p))/2, {l—p'{p))/2] can be regarded as "instanta- 
neous" critical points corresponding to the evolving crit- 
ical point (ii) introduced in Sec. II. In fact these points 
are on the line x + y = 1 and all trajectories approach 
this line as illustrated in Fig. 13 We note that along the 
line X + y = 1 the energy density of the barotropic fluid 
is vanishing {Qm = 0). 



C. p'(p)^-l 

Finally we consider the situation in which p'{p) ap- 
proaches — 1 asymptotically. One example is provided by 
the model (|37|l with 3/4 < < 1, in which case there is a 
type I singularity. Since p is given by Eq. H45|) in the limit 
oo, wehavep'(p) - -l-A(2/3-l)p2(/3-i) ^ -i-O 
and w 1 - y4p2(/3-i) ^ _i ^ o for A > 0. 

We study the case in which both A and B are positive. 
When x < the allowed range of x and y is the same as 
in the case of the subsection B. We show in Fig. ^the 
trajectories of the solutions for (3 — 0.85, A = 2 and B — 
1. The solutions eventually approach the point {x,y) = 
(0, 1) along with the line x + y = 1. This indicates that 
the trajectories evolve on instantaneous critical points 
{x,y) = [(1 +p'{p))/2, (1 -p'(p))/2] with the increase of 
p'{p) toward p'{p) — —1. The final state of the universe 
is dominated by a cosmological constant with vanishing 




-3 -2.5 -2 -1.5 -1 -0.5 

X 

FIG. 3: The phase plane for the model ^ with /3 = 1.1, 
a = 2f3 ~ 1, A = 2 and B = 1. The solutions approach 
"instantaneous" critical points: x = {1 + p' {p))/2 and y = 
(1 — p'(p))/2, which diverge as a; — > — oo and y —> oo ai the 
type III singularity. 



2.0 



0.50 



0.0 




-1.2 -1 -0.8 -0.6 -0.4 -0.2 

X 



FIG. 4: The phase plane for the model ^ with /3 = 0.85, 
a = 2(3 — 1, A — 2 and 5 = 1. The solutions approach the 
type I singularity at which the equation of state is u) = — 1 
with a; = and y = 1. 
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VI. COUPLED PHANTOM SCENARIOS 



The following EOS for dark energy is assumed: 



In this section we shall investigate a situation in which 
dark energy is coupled to dark matter. The coupled 
quintessence scenario was originally proposed in Ref. 1161 
as an extension of nonminimal coupling theories |l7j . 
This scenario has an advantage to explain the coincidence 
problem, that is, why the energy density of dark energy 
is comparable with the energy density of dark matter or 
usual matter. In Ref. jl^ the coupling between a scalar 
field (j) responsible for dark energy and ordinary matter 
is assumed to take the forms; Tj^^^^y^ ~ CTi.r„^(i)-u and 



T'', , = -CT,r 



where T, 



and T, 



are the 



energy momentum tensors of 4> and ordinary matter, re- 
spectively. Recently, a system of a phantom scalar field 
coupled to dark matter was studied in Ref. 18] by writ- 
ing the FLRW equations in an autonomous form, but 
this is not equivalent to our analysis since we assumed 
the EOS of dark energy rather than the existence of a 
scalar field. Furthermore the forms of the coupling Q are 
different from the ones we consider in this paper. 



A. Analytic solutions for coupled phantom 

Let the dark energy be a fluid whose equation of mo- 
tion is given by Eq. ^ . In general the coupling Q can be 
the function of a, H, p„n Pm, P and p. We can express 
the coupling Q that corresponds to scaling solutions, that 
is, the ratio r = Pm/ P is a constant. By Eqs. Q and Q 
we find 



Q , Q 



H 3H{wm — w) 



(64) 



We obtain the coupling Q for the existence of scaling 
solutions 0,123, as 



w 



PPr, 



Pn 



(65) 



which is derived by setting r — 0. In this case the energy 
density of ordinary matter does not vanish asymptoti- 
cally. 

Several authors considered the coupling of the forms 
like Q — 3HS{p+pm) |20] and Q = SHp,n 21], where 6 is 
a coupling strength. In this work we are interested in the 
coupling with which singularities discussed in preceding 
sections can be obtained analytically. 

Let us consider the following coupling as a toy model: 



(66) 



where (5 is a constant. As usually dark matter is regarded 
as a dust, i.e.. 



Pm = 0. 



(67) 



(68) 



that is, f{p) = pin Eq. (PT|) . 

Then combining Eq. |@J with Eqs. ^ and (01, we find 
a solution 



H 



2 /I 1 



3 V ^ ts-t 

\ t ^ ts - t ) t ' 
4/1 1 \ 1 



where 
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(69) 
(70) 
(71) 

(72) 



The same solution can be obtained if, instead of Eq. H66|l . 
the following coupling is adopted 



2 (p -f- Pm) 



(73) 



In this case appears as an undetermined constant. The 
solution (|69|l is the same as Eq. (|26|l with the choice n = 
2/3. Then the above solution gives rise to the Big Rip 
singularity at t = is. 

In the absence of the coupling Q, it is rather difficult to 
realize a model in which the equation of state w changes 
from w > —1 to w < —1 unless the function f{p) is 
double-valued as discussed in Sec. III. However this can 
be realized if the phantom is coupled to dark matter with 
an appropriate coupling. We note that the ratio r = 
Pm/p behaves as r = {tg ~ t)/t by Eqs. (|7n|l and itTTIl . 
which means that this ratio is dynamically changing and 
vanishes at t = is. Therefore the solution H69|) does not 
correspond to scaling solutions. This is understandable, 
since the system is completely dominated by the energy 
density of dark energy around the Big Rip. 



B. Multiple scalar field model 

It is possible to construct a scalar-field model which 
reproduces the solutions (|69() - (|71|l . This can be realized 
by a four-dimensional Lagrangian with two scalar fields 
<f> and X'- 

>C = ^i? - i(V0)^ + \{Wx? F(0, X) , (74) 
where the effective potential is 



if -<j>l -l± ^ rf +4>l 



ZtQ 



2f2 



Here (j)o and to are constants, and 



2 "^2 14 



e *o.(75) 



(76) 
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In Eq. (|74|l the sign of the kinetic term of x is oppo- 
site compared to the canonical one, which tells that the 
field X is regarded as a phantom. Then in this model, 
if we consider the quantization, the problem of nega- 
tive norm states would occur as the phantom is ghost 
miHil' When we adopt the Brans-Dicke type the- 
ory as in Ref . f^l , it is possible to avoid the problem of the 
negative sign in the kinetic term. Since it makes, how- 
ever, all things complicated, we pay our attention to the 
scenario given by 174|l in order to show how the solutions 
(|69|I - H71|I are realized in a simple scalar field model. 

The total energy density pt corresponding to the La- 
grangian (|74ll is given by 



Pt 



2t^ 



iie *o + ' 



2t^ 



-e 



(77) 



In the FRW background the equations of motion for 
scalar fields are 



-I- iH(t) - 
X + SHx 



Tl _±±x 

-e 't'o 



n2 



■*o -I- . ' 



'-e-^ =0, (78) 

00*0 

O-e ^0. (79) 



Then by using the Friedmann equation — we 
find a solution: 



H = 



2 [ t 

t 



1 



t.-t 



</) = (/)o In -- , X 00 In 

to to 



(80) 
(81) 



where tg appears as a constant of the integration. Then 
if one chooses 00 as 



p as 



1 :^ 



e *o 



2t^ 



2t^ 



'2^ 2tf 











2tl 



(84) 



it follows 



2,H {prn + Pm) = - {p + 3H {p + p)} 



2(j)oti 



4-+X 

-e 



(85) 



By comparing Eq. H85|l with Eqs. Q and one gets 



Q 



24>otl 



+ X 



(86) 



Note that all the parameters in the action (|74|l have 
the scale of order the Planck mass. The parameter ts ap- 
pears dynamically as a constant of the integration. If the 
present universe corresponds to the time t of order ts , one 
can estimate as ts ~ 1.8 x 10^*^ years ~ (lO~^'^eV) ^ . 
The coupling Q in Eq. H86|l cannot be straightly rewritten 
in the form of Eq. or Eq. but with the solu- 

tion, all the expressions for H , pm and p give the same 
functions in terms of the time t. 

One can consider a generalization of the model in 
Eq. |(73I) as follows: 



1 
2^2 



1 



1 



(V0)2 + -(Vx) 



ri9 



2f2 



■24> 
<^0 



xl 



2tl 



(87) 



Here 0Oj XOj and to are constants and 



4 
3^ 



(82) 



^ 2 i4 

2'* -^0 > 



3 2 4 



the Hubble rate H in Eq. H69|l is reproduced. 

The field is identified as dark matter with an energy 
density pm and the field x as dark energy with an energy 
density p. The total energy density pt in Eq. (|77|l may 
be given by 



1 T -^ttX 

-o -I ^e '*o 

2^ 2tl 



2tl 



i>o 
-e *o 



^/l 
2 V* 
1 



1 



ts 



t / t 



—X H — ''■o -t- 

2 



2 



1 



1 



t ts-t t 



1 



(83) 



which agree with Eqs. H7U|) and (|71|l by using Eqs. H7t)|) 
and (|82|l . Furthermore if we define the pressures pm and 



Then a solution has the following form: 



2 V t 



xl 

ts-t 



= 00 In — , X = xo In — — 

to to 



(89) 
(90) 



The Hubble parameter has a minimum a,tt — 0ots/(0o + 
Xo) for Xo/00 > and diverges for i — > and t ts- 
When t <C is the system is dominated by the field 
and the effective EOS is given by w ~ — 1 + 4/(3k20q). 
On the other hand the field x dominates around t ^ 
ts, corresponding to the equation of state: w ~ — 1 — 
4/(3k2xo). Then the system transits from the region 
w > — 1 to the region w < — 1 and approaches the Big 
Rip singularity at t = is. 

We may investigate whether the present model can sat- 
isfy the present observational data. The deceleration of 
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the universe changed to the acceleration about five bil- 
Uon years ago. Let us write the corresponding time as 
t ~ tc- Since 

a/a = H^ + H 

4 I t2 

t[ts-t) / ' 

then tr satisfies 



{ts - ty 



(91) 



{c^l~xl){^'{<l>l-xl)-2]tl (92) 
-2^1 {k^ {4>l - xl) - 2} + {K^cj^l - 2) cgtl = . 

According to the present data, the density of dark energy 
is 72% of the total density and that of dark matter is 21% 
[25| . We define the energy density of dark matter pm and 
that of the phantom p by 



2^ 



2 , Tie 



e *t> H — ^ 



e *o 



P = —^x 



t ts-t 

riO x_ 

—g-e *o xo H — ^ 



+ XO ^-2 



e 



xl 



A{ts~t)\t ts-t 



(93) 



If t corresponds to the present age of the universe, that 
is, fourteen billion years, one finds 



tx 



{ts - t) ^0 



Then 



2x| 
70g 



1 



tc 



14 ^-irl^ + i"* 



14-5^ 9 ^2xo 



14 V702 



If we define a coincidence time t by p = p™, one has 



(94) 

(95) 
(96) 

s 

(97) 



Dark matter is regarded as a nonrelativistic matter 
with EOS: Wm = 0. We may consider a situation in 
which this originates from the field (p with constant EOS: 
m„i = -1 + 4/(3k2(/)2), which gives = 4/(3^^). We 
caution that this property does not necessarily hold, since 
this constant EOS does not correspond to the attractor 
solution as shown in Sec. V. We just wish to estimate the 
approximate time tg by imposing this condition. Then 



substituting 
obtain 



4/(3k2) and Eq. ^ for Eq. we 



K^xl = 0.0831097., 



(98) 



This gives ts — t ^ Q.02t ~ 0.3 billion years toward the 
Big Rip, which is rather a short time. The coincidence 
time t is given by t — t ^ 0.04i ^0.6 billion years ago. 



VII. ACCOUNT OF QUANTUM EFFECTS 
AROUND THE SINGULARITIES 

It is quite natural to investigate the role of quantum ef- 
fects near singularities. Indeed, the dark energy universe 
with a singularity typically evolves to its end with the 
growth of the energy. As a consequence the curvature of 
the universe grows as well, which implies the beginning 
of a second quantum gravity era. The type IV singular- 
ity introduced in Sec. Ill is not of this sort because the 
curvature does not grow significantly there. 

In what follows we shall study the effect of quantum 
backreaction of conformal matter around the type I (Big 
Rip), type II (sudden future) and type III singularities. 
In these cases, the curvature of the universe becomes 
large around the singularity at i = t^, although the 
scale factor a is finite for type II and III singularities. 
Since quantum corrections usually contain the powers 
of the curvature or higher derivative terms, such correc- 
tion terms play important roles near the singularity. We 
now include the quantum effects by taking into account 
the contribution of conformal anomaly as a backreaction. 
The conformal anomaly Ta has the following form: 



TA^blF 



-DR 
3 



b'G + b"UR 



(99) 



where F is the square of a 4d Weyl tensor, G is a Gauss- 
Bonnet curvature invariant, which are given as 



F = (l/3)i?2 - 2i?,jff^ 



(100) 



In general, with N scalar, N1/2 spinor, A^i vector fields, 
N2 (= or 1) gravitons and A^hd higher derivative con- 
formal scalars, the coefficients h and b' are given by 

_N + 6A^i/2 + l2Ni + &IIN2 - 8AfHD 
120(47r)2 ' 

A + llA^i/2 + 62jVi + 1411A^2-28A^HD , . 

360(47r)2 ' 

We have 6 > and 6' < for the usual matter except 
for higher derivative conformal scalars. Notice that b" 
can be shifted by a finite renormalization of the local 
counterterm R^ ^ so b" can be arbitrary. 

In terms of the corresponding energy density pA and 
the pressure density pa, Ta is given hy Ta = — Pa + ^Pa- 
Using the energy conservation law in FLRW universe: 



PA + 3H {pA +Pa) = 0, 



(102) 
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we may delete pa as 

Ta = -4:PA - Pa/H . 
This gives the foUowing expression for pA'- 

PA ^ j dta^HTA 



(103) 



dt a^H 



UbH^ + 2Ab'{-H^ + H^H + H 



{4b + 6b")(H + 7HH + 4:H'^ + 12H^H] . (104) 



A different form of pA was obtained in Ref. by re- 
quiring that the quantum corrected energy momentum 
tensor Ta^^ has the form as TAfu, = {TA/'^)g^iu in the 
conformal metric case rather than assuming the conserva- 
tion law H102|l . The quantum corrected FLRW equation 
is 



—H^ = p + PA- 



(105) 



Since the curvature is large around t = tg, one may as- 
sume (3/K2)iJ2 < \pa\- Then p ~ -pA from Eq. (fTnijl . 
which gives 



p + AHp 



= H 



-l2bH^ + 246'(-ij^ + H^H + ij") 



-(46 + 66") ( H + IHH + AH"^ + UH'^H 



(106) 



From the energy conservation law p + 3H {p + p) ~ for 
p = — p — f{p), one obtains 



H 



P 



3f{p) 



(107) 



A. Type III singularity 

For the type III singularity, the energy density p di- 
verges att = tg. Let us consider the function f{p) ^ Bp^ 
discussed in Sec. IV. In the absence of quantum correc- 
tions the type III singularity appears for (3 > 1 and p 

2 

behaves as p oc (ts — t) ^-"^'^ . When quantum corrections 
are taken into account, it is natural to assume that near 
the singularity p behaves as 



p = Pq (is - 



(108) 



As p may diverge at i = ig, we consider negative values 
of 7. Using Eq. I|107(l . one finds 



- 1-/3 

3B 



(109) 



Since we are considering the case /3 > 1 and 7 < 0, 
we have that 7 (1 — /3) > 0. By picking up most singular 
term in the r.h.s of Eq. H106|) . it follows 



-6 ( |6 + 6" ) UR. 



(110) 



Then substituting Eqs. lfTU5|) and ifTU^ for Eq. ((TTn|l . we 

obtain 



7 



1 - 2^ 

This means that p and H evolve as 



p(x{ts- t)i-2 



(111) 



(112) 



around t = tg. Compared with the classical evolu- 

2 

tion: p (X {tg — t) , the energy density diverges 
more rapidly because of the condition: 4/(1 — 2/3) < 
2/(1 - 2/3) < for /3 > 1. On the other hand the Hub- 
ble rate behaves as _ff cx {tg — t) 1-2/3 jn the classical 
case, which means that H is less singular in the pres- 
ence of quantum corrections because of the condition: 
(3 - 2/3)/(l - 2/3) > 1/(1 - 2/3) for /3 > 1. 

We numerically solve the background equations and 
show the evolution of the Hubble rate for /3 = 2 and 
_B > in Fig. |S1 In the presence of quantum corrections 
one has H (x {tg — ty^^ around t = tg^ which means 
that H approaches zero. Meanwhile in the absence of 
quantum corrections we have H oc (tg — t)~^/^, thereby 
showing the divergence of H at t — tg. These properties 
are clearly seen in Fig. [3 

From Eq. Hl()9|l we obtain 



a ~ oo exp 



1-/9 

Po 



35(1-/3) 



{tg - t) 



7(1-/9) 



(113) 



where oq is a constant. Comparing the classical case [7 = 
2/(1 — 2/3)] with the quantum corrected one [7 = 4/(1 — 
2/3)], we find that the power of {tg — t) is larger in the 
presence of quantum corrections. Then the scale factor 
approaches a constant ao more rapidly if we account for 
the quantum effect, implying that the spacetime tends 
to be smooth, although the divergence of p is stronger. 
Thus quantum effects moderate the classical singularity. 



B. Type I singularity 

We shall next consider the model f{p) ^ Bp^ with 
1/2 < /3 < 1 when p is large. In this case there exists 

the Big Rip singularity as shown in Sec. IV. We note that 

2 

the classical evolution is characterized by p cx (tg — t) 1-2/3 

and H (X {tg ~ t) , both of which exhibit divergence 
for /3 > 1/2. 

When the quantum correction is present, let us first as- 
sume that the time-dependence of p is given by Eq. H108|) 
with negative 7. Since 7(1-/3) < in this case, we 
might expect that Eq. H106|) would give the following ap- 
proximate relation around t — tg: 



6b'H^ 



(114) 



The term on the r.h.s. grows as iJ* cx {tg — i)^4+47(i-/3)^ 
but this does not give a consistent result, since p becomes 
negative for 6' < 0. 
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(b) 


/' 







0.5 



t 



1.5 



FIG. 5: The evolution of H for the model /(p) = Bp^ with 
l3 — 2 and B > 0. The case (a) corresponds to the one 
in which quantum corrections are taken into account with 
coefficients b = 0.5, b' = —0.1, and b" = 0, whereas the 
case (b) does not implement such effects. The Hubble rate 
approaches H = Q with a finite time in the case (a), while it 
diverges in the case (b). 



This tells that our assumptions should be wrong and p 
does not become infinite. If p has an extremum, Eq. H1U7|) 
tells that H vanishes there since p — Q. In order to 
confirm whether the quantum effect moderates the sin- 
gularity or not, we have numerically solved the back- 
ground equations for the model where p is exactly given 
by f{p) = Bp^ with 1/2 < /? < 1 and coefficients 6 > 0, 
6' < and b" = 0. We find that the Hubble rate ap- 
proaches zero with a finite time, as is similar to the case 
(a) in Fig. |S1 Thus the presence of the quantum correc- 
tion moderates the Big Rip singularity as well. 



C. Type II singularity 

Finally we consider the model f{p) = C (po — p) ^ 
with 7 > as an example of the sudden future type sin- 
gularity. By Eq. ^ the energy density p in the classical 
case behaves as 

p^po- {kCV3a^(7 + 1) {ts - t)y*' , (115) 

around the singularity ai t — tg. Here we write to in 
Eq. (|62l) as tg. Using the Friedmann equation = 
we find 



H 



1 



kC./3^{j + 1) its - t) 



Since < 1/(7 + 1) < 1, diverges at t = ts while H is 
finite there. 

Let us now include the quantum corrections. We as- 
sume the following form of p around t ^ tg'. 

p^ pa + pi{ts-t)\ (117) 

with a positive constant v. Using Eq. H107() . one gets 



H 



1+7 
3{-C) 



(tg-t) 



-H- 1/(1+7) 



(118) 



Since (1 -f 7) > 0, we can use the same approximate 
equation as in Eq. H11U|I . This gives 



27 + 1 ' 



and 



P = Po + Piit 
H cx {ts 



t)2T+l 



(119) 

(120) 
(121) 



(116) 



Since (274-3)/(27-|- 1) is larger than 1, not only H but 
H are finite in the presence of quantum corrections. This 
is numerically confirmed in Fig. Thus it is clear that 
quantum effects work to make the universe less singular 
or completely non-singular (basically, asymptotically de- 
Sitter). It was shown in Ref. ^ that this property also 
holds for scalar-field dark energy models. 

Finally we should mention one point. When the quan- 
tum correction becomes important, this typically works 
to provide a negative energy density pA which nearly can- 
cels with the energy density p of dark energy. This is the 
reason why the Hubble rate does not diverge if we ac- 
count for the quantum effect. 



VIII. DISCUSSION 

In summary, we discussed the fate of the (phantom) 
dark energy universe by assuming the equation of state 
(EOS) with a form: p = —p — f{p). Our main interest 
is to clarify the structures of future singularities which 
appear with finite time ts ■ We classified the types of sin- 
gularities into four classes. The type I is the Big Rip 
singularity at which all oi a, p and p exhibit diver- 
gences. The type II corresponds to the sudden future 
singularity in Q at which a and p are finite but p di- 
verges. The type HI is similar to the Big Rip, but the 
scale factor is finite. The type IV is a mild singularity 
at which a, p, and p are finite but higher derivatives of 
the Hubble rate diverge. The model given in Eq. H32(l 
includes all of these singularities, which are investigated 
in details in Sec. III. In Sec. IV we adopted simpler forms 
of /(p) which are obtained by considering limiting cases 
of the function (|32|) and studied the general structure of 
singularities. This is very helpful to understand when 
these finite-time singularities appear. 
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FIG. 6: The evolution of H for the model /(p) = C(p-po)"^ 
with 7=1/2 and C < 0. We implement quantum effects in 
the case (a) with coefficients 6 = 0.5, b' — —0.1 and b" = 0, 
whereas the case (b) does not account for such effects. We 
find that H approaches zero with a finite time in the case (a), 
while it diverges in the case (b). 



In Sec. II we wrote the background equations in an 
autonomous form and showed that there exists one sta- 
ble critical point if the constant EOS {w = p/p) of dark 
energy is less than — 1 . It was also found that there is no 
stable node for w > —1. In Sec. V we have numerically 
confirmed this property for the function (|30|l that asymp- 
totically gives constant values of w. The singularity for 
w < — 1 with the function H3(J|I corresponds to the type 
I, which implies that the Big Rip singularity is a late- 
time stable attractor. We also performed phase-space 
analysis for the model (|32|l and found that the solutions 
approach "instantaneous" critical points introduced in 
Sec. II. We have numerically checked that the late-time 



solutions agree well with analytic estimations. 

In Sec. VI a (phantom) dark energy scenario coupled 
to dark matter is investigated. We constructed the form 
of the coupling Q that reproduces the background evo- 
lution (|26|l giving the Big Rip singularity a,t t = tg, see 
Eqs. and H73|l . Unlike the case of a single fluid, 
this allows a possibility to lead to the transition from 
w>— ltow<— 1 without using a double- valued func- 
tion for f{p). We also explicitly presented a two-scalar 
field model which gives the same singular behavior for 
the background evolution. 

Finally, we accounted for the quantum correction com- 
ing from conformal anomaly and studied its effect to type 
1, II and III singularities. We found that this typically 
works to moderate the singularities or prevent them by 
providing a negative energy density. In this case the di- 
vergent behavior of the Hubble rate is prevented by the 
quantum effect, which is followed by the decrease of H 
toward 0. Therefore the presence of the quantum effect 
can drastically change the evolution of the universe. For 
instance, the new inflationary era in far future may be 
possible according to the conjecture of Refs. |E Sl- 
it is quite clear that with the growth of the energy 
density near to singularity the effects of quantum gravity 
(string/M-theory) may become dominant. Hence, it will 
be extremely important to estimate the effect of dilatonic 
or modulus higher-order corrections 27] in dark energy 
models motivated by string theory, which we leave to 
future work. From another point, even if stringy effects 
stop the future singularity and possible phantom phase 
is transient, the growth of dark energy density may give 
its imprints to current universe. It would be a challenge 
to search for such imprints. 
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